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Abstract 

The computer language Logo facilitates the teaching of analytic geometry and calculus from the notion of 
curvature, through its "turtle geometry" facility [2]. We provide some Uicorctical basis for finding turtle 
geometry equivalents of familiar curves in analytic geometry, and vice versa, by some simple mcUiods 
apparently previously unnoticed. In particular, we study turtle geometry programs where the curvature of a 
line is a trigonometric flmction of its orientation. 

1 . Introduction 

Recently the computer language Logo has become available for the first time on a wide variety of personal 
computers (1, 6, 3). A distinctive feature of Logo is an emphasis in the graphics package on "turtle geometry", 
a geometry emphasizing discrete changes to local curvature of curves rather titan coordinate systems. Despite 
its simplicity, turtle geometry is quite powerful, and provides an alternative viewpoint that makes certain 
mathematical concepts clearer, as discussed in [2). 

But an issue not well addressed in [2] is relationship of curves created by turtle geometry to the familiar 
ones of analytic geometry. In tliis paper, thus, we present metliods for finding small-step turtle-geometry 
equivalents (mostly previously unnoticed) for some well-known analytic-geometry curves, and vice versa. 
This material provides a good exercise for calculus students, while also demonstrating some numerical 
analysis concepts. 

2. Turtle geometry 

Turtle geometry draws planar curves using two commands, FORWARD and RICH L A cursor called a 
turtle represents a point and orientation on a graphics screen. FORWARD causes the turtle to move forward 
in a straight line a distance (approximately in millimeter units) specified as its argument, in die direction the 
turtle is pointing. RIGHT causes the turtle to turn in place a number of degrees specified by its argument. 

Turtles tluis have a state consisting of a screen position and an orientation (or heading). The Cartesian 
coordinates are referred to as XCOR and YCOR, and the heading as IIHADING (though we will use the 
abbreviation H for it in this paper). Wo assume the turtle starts with XCOR = 0, YCOR = 0 (the center of the 
screen), and FIFAD1NG = 0 (pointing straight up). And a IIHADING of 90 is straight to die right, 180 
straight down, and 270 straight to the left. Note that 1 IHADING can be thought of as the toul amount a 
turtle has turned right since it started. 

Logo piograms are written in terms of procetlures. A procedure consists of a sequence of commands, on 
successive lines. Ihe first line contains the word TO followed by die name of the procedure, and then 
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arguments if any. 1 lie last line of a procedure is signalled by I-ND. 

3. Curvature 

We start from the classical definition of curvature in a Cartesian coordinate system, as given in most 
advanced calculus treatments (e.g. [4], p. 305), assuming that clockwise is positive curvature: 

K(x,y) = -(d^y/dx^)/ [1 + (dy/dx)^]' ^ 

We can put tliis in a fonn closer to turtle geometry by noting that dy/dx = cot(H), H the heading of the turtle 
[5]: 

K(x,y,H) = -(d^y/dx^) / [I + cot^H]^ ^ = -(d^y/dx^) * sin^H 
We call tills the "hybrid" curvature formula because it includes both Cartesian and heading variables. 

Curvature is die rate at which a curve turns clockwise, and thus can be approximated in Logo by die ratio 

of the amount a turtle turns RIGHT to the amount it goes I'ORWARD -- provided the turtle takes small 

enough steps so as to make its padi look sufficiently smooth. (Practically speaking, this means something like 

steps of 1 millimeter and turtle turns of five degrees on most graphics terminals that use Logo.) So if we 

define the curvature function we want a curve to follow by die procedure K of no arguments, cidier of these 

recursive Logo procedures will probably draw what we want: 

TO CURVE 
FORWARD 1 
RIGHT K 
CURVE 
END 

TO CURVE 

FORWARD (QUOTIENf 1 K) 

RIGlif 1 

CURVE 

END 

We shall call diem the "angle-control" and "step-control" variants, respectively. Which should we use? 
Unfortunately, sometimes they do not work die same, and only one gives to the correct result. This is because 
these arc discrete approximations of smooth curves, and anytime the absolute value of K becomes large in the 
first version, or the absolute value of K approaches 0 in the second, the approximation bicaks down. So we 
will have to watch for that. 

In what follows we shall only concern oui'sclves with the shape of a curve, and not its si/c or orientation. 
I bus we shall ignore constant muliiplicrs in front of a curvature formula, except sometimes the sign, and we 
shall throw away any constant midtiplicrs generated in analysis routinely. We ;ilso will not concern ourselves 
with at what screen k)cation we start drawing the sh.ipe (e.vccpt in section 4.3), since Un tie geometry docs not 
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use a coordinate system, though the usual convention in Logo is to start at tlie center of a grapliics screen. But 
it will make a difference what heading we start with in the procedures we discuss -- different pieces of the 
curve will appear for different starting headings, and for some headings not defined for the curve (e.g. 
asymptotes), pathological curves like straight lines will be generated. 

Special note: in this paper we assume that headings of greater than 360 degrees or less than 0 are 
meaningful and significant as indications of "how far the turtle has turned". Most Logo implementations, 
however, take "heading" as this modulo 360. As a result, some of the Logo procedures discussed below will 
not execute properly in those implementations (in particular, those where heading (H) appears other than as 
an argument to a m'gonometric function). But it is easy to define a counter variable to use instead of 
HEADING. 

4. Some examples: from y = f(x) to K = g(H) 

We give now some examples of curves initially described as y = ffx) in new terms of their curvature as 
function of heading, K = g(H). The method is (1) find dy/dx and d^y/dx^ in terms of x; (2) express d^y/dx^ 
in terms of dy/dx only, no x’s; and (3) then substitute cot(H) for dy/dx, and multiply by -sin\H) to get the 
curvature K in terms of H. I'his metliod is not guaranteed to work because step (2) is not possible for some 
f(x), but it will work for some interesting functions. 

4.1 . Conic sections 

Let us apply tlie hybrid foimula to some familiar curves. First, the general parabola y = ax^ + bx + 

c. Since d^y/dx^ = 2a, a constant, following our policy of ignoring constants we just turn right .siii^(H) at each 

turn in the "angle-control" procedure, as follows: 

TO PARABOLA 
FORWARD I 

RIGH r (PRODUCT (SIN HEADING) (SIN HEADING) (SIN HEADING)) 

PARABOLA 

END 

If the turtle is started with a heading between 0 and ISO, it will draw a downward-tending parabola; else, an 
upwards-tending one. 

One of the first things people learn to do with Logo is to draw a circle by repeating many times a 

forward-right pair of commands with constant arguments: 

TO CIRCLE 
FORWARD 1 
RIGHT 1 
CIRCLE 
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END 

( I'liis is a special ease of the well-known Logo procedure, F’OLY.) We can prove tliat in fact y = -/(r^-x^) is 
drawn with this procedure, dy/dx = -x/7(r^-x^), and d^y/dx^ = r^/(7-x^)^^ Hence 1 + (dy/dx)^ = 
r^/(r^-x^), and (1 -f (dy/dx)^)’ Vr = d^y/dx^. Substituting cot(H) = dy/dx, d^y/dx^' = (1 -f cot^(H))’^'^ = 
csc\H). Hence die curvature is -sin^(H) * csc^(H) = -1, a constant. 

We can generalize diis to ellipses and hyperbolas. Hirst, consider die general ellipse y = cV(a^-x^), where c 

is die ratio of die height to the the width. Then dy/dx = -ex/ 7(a^-x^) and d^y/dx^ = caV(a^-x^)'^ ^ Hence 

1 + (dy/dx)Vc^ = a7(a^-x\ and c(l -f (dy/dx)Vc^)^'^ = d^y/dx^ = c(l -f (l/c^)cot^(H))*-^ = (l/c^)(c^ 

+ cot^(H))^'^ Multiplying by -sin VO. and ignoring die multiplicative constant, diis then leads to a total 

curvature of -{c^sin^(H) -f- cos^(H))^'^ = -[(c^-l)siii^(H) -f- 1]’ ^ Since sin^(H)>0 and c^-l>-l forc5*0, the 

curvature is always negative for cv^O, but the degree of curve varies periodically with the heading. Here is the 

l.ogo procedure for diis general ellipse: 

■fO RLLIPSF. :C 
FORWARD 1 

RIGH'f THRFFHALVES DIFFFRENCF: (PRODUCT (DIFFERENCE SQUARE :C 1) 

SQUARE SIN HEADING)! 

ELLIPSE :C 
END 

TO THRFiEHALVES :X 
SQR r (PRODUCT ;X :X :X) 

END 

By similar analysis we can get the curvature of a general hyperbola, y = c7(a^ + x^). It turns out dy/dx = 
cx/7(a’-l-x^) and d^y/dx^ = caV(a^-l-x^)^^ so c(l - (F/c^ot^(H))''^ = d^y/dx^ and the curvature is 
proportional to -(c^siir(H) - cos^(ll))^^ = -((c^-F l)sin^(H) - 1]''^. Ihis has maximum curvature at H = 90 
degrees, and dien gradually turns left more and more slowly, reaching an asymptote when the curvature 
becomes zero, or (c“+ 1 )sin^(H) = F, or H = arcsin[ 7(I/(I -i-c^)] = arccot(c). 

Note both the ellipse and hyperbola become a parabola when c = 0, since (l-sin^H))’ ^ = cos^(H). And 
when c approaches oo, both approach (c“sit7(l I))' ^ = c\in\H), which is also a parabola. 

4.2. Inverses 

Take f(x)=e’‘; dy/dx = d^y/dx^ = e*, so d^y/dx^ - cot(H), and the curvature needed to draw it is 
-sin^(H) * cot(H) = -sii7(H)cos(H) = -.5 sin(l!) sin(2H). 



Since the inverse function of f(x) looks just like it but "turned over" and oriented in a different direction. 
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we should expect their curvature formulas to be related. We can verify this for the inverse function of c’', 
ln(x); dy/dx = 1/x and d^y/dx^ = -1/x^, so d^y/dx^ = -cot^(H), and the curvature is sin(H)cos^(H) = .5 
cos(H)sin(2H). 

In general, we can obtain the curvature of the inverse of a fiinction by interchanging all sin(H) and cos(H) 
terms and then reversing the sign of the total expression so that it curves in the opposite direction. This is 
because 

dx/dy = 1 / (dy/dx) = tan(H) 

So tan(H) should be substituted forcot(H) in all occurrences of it, and vice versa. But since 

1 + (dy/dx)^ = 1 + cot^(H) = l/sin^(H) 

1 + (dx/dy)2 = 1 + tan^(H) = l/cos^(H) 

sin(H) and cos(H) must necessarily interchange too. 

4.3. Problems with inflection points 

Unfortunately, we usually cannot use an expression for die curvature in terms of only die heading when die 
curvature changes sign somewhere. Often diese conditions arise with square roots which can be cither 
positive or negative. An example is y = x^ dy/dx = 3x^ and d^y/dx^ = 6x. So trying to express d^y/dx^ in 
terms of dy/dx, we get ± V[12(dy/d.x)] = ± V[12cot(H)]. and a curvature of Tsin^(H) V[12cot(H)]. 

But we cannot choose only one sign to cover die whole range - die sign* must change at the inflection point, 

the point where 6x = 0, x = 0, or tlie place where arccot(3(0) ) = H = 90. We must tlicrefore add a 

conditional branch, something like tliis (assuming we start with a heading which is a small positive number); 

TO CUBE 
FORWARD 1 

RIGITI PRODUCT (SQUARE SIN 11 HADING) (SQRT COT HEADING) 

IF HEADING = 90 [CUBERIGHTSIDE] ELSE [CUBE] 

END 

TO CUBERIGHTSIDE 
FORWARD 1 

RIGl rr MINUS PRODUCT (SQUARE SIN HEADING) (SQRT COT HEADING) 

CUBI-RIGHISIDE 

END 

where SQR T is the square root function (sec the Logo manuals cited for the exact interpretation of this 

syntax). But since tlie curve is really composed of small straight lines, with some roundoff error in calculating 

the heading over a period of time, we may never actually reach a heading of zero. Instead we had better say 

for procedure CUBE: 

TO CUBE 
FORWARD I 
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RIGHT PRODUC T (SQUARH SIN HHADING)(S0RTCO T HHADING) 
IF HHADING < 89 [CUBE] EI.SE(CUBERIGH TSIDE) 

END 

with CUBERIGH'l'SIDE remaining as above. 



5. Parametric equations: from x = f^(t), y = to K = g(H) 

Some curves arc easier to express in terms of parametric equations for x and y tlian in tenos of a equation 
including both x and y. To find curvature in terms of heading for them tlicre arc two steps: (1) express the 
parametric equations in tenos of H entirely (use dy/dx = cot(H)), and (2) compute derivatives with respect to 
tliis expression to obtain the curvature. We can rewrite die basic curvature fonoula in terms of derivatives of 
H: 

K(H) = -sin^(H)*dVdx^ 

But d^y/dx^ = d/dx (dy/dx) = d/dx (cot(H)) = csc^(H) * dH/dx 
Hence K(H) = -sin(H) / dx/dH 

But also 

dH/dx = dH/dy * dy/dx = dH/dy * cot(H) 
so also K(H) = -cos(H) / dy/dH 



As an example, take the cycloid described by 
X = a(t - sin(t)), y = a(l - cos(t)) 

cot(H) = dy/dx = (dy/dH)/(dx/dH) = sin(t)/(l -cos(t)) = cot(t/2) 

and hence t/2 = H, t = 2H, and 

X = a(2H - sin(2H)), y = a(l - cos(2H)) 
dx/dH = 2a(l -cos(2H)), dy/dH = 2a sin(2H) 

So tlic curvature is, using die simpler expression dy/dH: 

•sin(H) / 2a sin(2H) = -l/4acos(H) 

and hence the figure can be drawn by die Eogo procedure 

TO CYCEOID 

FORWARD SIN HEADING 

RIGH T I 

CYCLOID 

TNI) 

Note we use die "step-control" procedure format because scc(H) becomes infinite at ± 180 degrees. 



6. Polar coordinates to K = g(H) 

We can also dr<iw curves given in polar coordinates, r = f(^); 

X = rcos{^), y = rsin(0) 

dx/d0 = (dr/d^)cos(0) - r sin((/), dy/dO = (dr/d(?)siii(^) -F rcos(^) 
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cot(H) = dy/dx = (^y/3^?)/(9x/^^) = 

[(dr/d^)sin(^) + r cos(0)] / [(dr/d^)cos(^) - r sin(0)) 

= cos[0 - Arctan[(dr/d^)/r]] / -sin[0 - Arctan[(dr/d0)/r]] 
= -cot(^ - Arccot[r/(dr/d0)]] 

Hence 

H = -^ 4- Arccot[r/(dr/d^)], or cot(H + 6) ■= r/(dr/d^) 



The curvature of a curve expressed in polar coordinates can be found by manipulation of die x-y curvature 

formula, and is given in many books (e.g. [4], p. 306); 

K(0) = [r^ + 2(dr/di9)^ - rfd^r/di?^)] / [r^ + (dr/dl9)Y ^ 

We substitute into it die rcladonship just derived, i.e. (dr/d^)/ r = tan(H + 0)'. 

K«?,H) = [1 + 2tan^(H + ^)-(dVd0Vr]/r(l + tan^(H + 

= [2cos(l I + 6) - cos\H + 6) - [cos^(H + 0) * (d^r/d<?^)/r]] / r 
= [2cos(H + ^)]/r- [cos\h + ^) * [1 + (d\/d^Vr^]] 

For particular situations we may be able to express the above only in terms of H. For example, take the 

logarithmic spiral r = for which tan(H + 0) = (dr/d^)/r = = c. Hence: 

dr/d<? = ce'"'’, d\/dd'^ = cV'^ 

K(H,r) = 2cos(arctan(c))/r - cos\arctan(c)) * [1 + c^]/r 

Which just factors into a constant term divided by r = e*^^ = eC(arccot(c)-n) curvature is just a 

constant times e'*^^ = b ", b = e^ So it can be drawn by this procedure: 

TO LOGSPIRAL 

l-'ORWARD EXPONENT 1.01 H 

RIGHT 1 

LOGSPIRAL 

END 

where liXPONBNT represent die first argument taken to the power represented by die second argument. 



7. Changing-argument recursive-procedure forms 

There's a simpler way we can express LOGSPIRAL (given in chapter 2 of [2]), by giving it an argument: 

TO LOGSPIRAL -.SIDE 
F0RWARD:S1I)E 
RIGHT 1 

L0GSP1RAL(:S1DE* 1.01) 

END 

where :S1DF, is a parameter to the procedure. 



Many procedures diat draw smooth curves c;m be written in this form. Letting the airvature on two 
successive recursive calls be and K^, we try to find a simple expression relating to Kj, preferably 
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independent of heading. In die above, K2=Kj/1.01 which is independent of heading. Note if die curve is 
smooth enough, we can approximate by a Taylor series about Kp that is 
KjCH) = Kj(H) + (SK(H)/6H)AH 

so dicrc is always an "additive" approximation for any Logo step-control procedure: 

TO CURVH :CURVATURE 

FORWARD 1 

RIGHT .CURVATURF, 

CURVE SUM .CURVA fURE {DERIVE HEADING) 

END 

where DERIVE is a procedure diat computes the derivative of die curvature function. Rut since this is only 
an approximation it will (usually) gradually diverge from the true curve, though slower for some functions 
than others. 



8. From K = g(H) to y = f(x) 

Now we address the inverse problem, which is to determine the Cartesian description of a function given 
its curvature as a function of heading. Two approaches arc possible. 



8.1 . From g(H) to y = f(x) 

Let i. = dy/dx. I’hcn if the curvature is a hmetion of heading g(H): 

g(H) = -sin^H) * dz/dx 
-csc^{arccot(z)) * g(arccot(/.)) = d/,/dx 
dx = dz. / [-csc\arccol(z)) * g(arccot(z))] 

X = -/dz/(g(arccot(z))(l-f z.^)' 



In particular eases we may be able to solve this, then invert it and integrate with respect to x to get an f(x). 

For instance, consider die curve drawn by the procedure 

TO MYSTERYl 

FORWARD 1 

RIGHT SIN HF:AD1NG 

MYSTERYl 

END 

which draws a kind of "squashed parabola". Its curvature is E(1 1) = sin(H). Hence the equation is 

X = -/dz./(l 1-/^) = C - arctan(z.) 

Soz - tan(C-x) = dy/dx. 

Hence y = / tan(C-x)dx = ln|cos(C-x)| 

which is thus die mystery curve as a function of x. 



As another example, consider the similar procedure 
TO MYS rFRY2 
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FORWARD 1 

RIGHT PRODUCT SIN HEADING SIN HEADING 

MYSTERY2 

END 

which draws a less-squashed parabola. Its curvature is K(H) = sin^(H), and: 

X = - f dz/ V([l + = C-arcsinh(z) 

Soz = sinh(C-x) = dy/dx. 

Hence y = / sinh(C-x)dx = cosh(C-x) 

which is a catenary. 



This metliod depends however on tlircc steps, two integrals and a function inversion. If any of these are 

not possible, the method will not work. But anytime the first integral is possible, one does derive a differential 

equation in tenns of dy/dx and x, and it may be possible to solve this with other methods, or perhaps 

approximate a solution by infinite series. P'or instance, tlie procedure 

TO MYSTERY3 
FORWARD 1 

RIGHT EXPONENT (SIN HEADING) 5 

MYSTERY3 

END 

leads to the differential equation 
(dy/dx)V3 + (dy/dx) = x 
which has no easy solution. 



8.2. From g(H) to x = f(H) and y = h(H) 

If, however, it is acceptable to determine llie curve in terms of the parameter 11, tlien we can often solve 
other situations too. I^ecall die formulas from section 5 for curvature in terms of tlie lieading parameter 
exclusively: 

K(H) = -sin(H) / (dx/dll) , or K(H) = -cos(H) / (dy/dH) 

We can rearrange tliese expressions for K(H) and integrate tliem: 

dx = -sindDdI I / K(H). dy = -cos(H)dH / K(H) 
x = -/[sinODdll / K(H)] 
y = -/[cos(H)dll / K(H)] 

Sometimes we can then find a substitution to eliminate H and get an equation in terms of only x and y, but 
not always. 

As an example, consider curves drawn by tlie "step-control" form 
TO AMPOLY:C 

1-OUWARD SIN PRODUCr :C HI-ADING 
RIGHT 1 
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AMPOLY ;C 
END 

where :C denotes a parameter. (AMPOEY stands for "amplitude modulation POLY", since the curvature of 

a circle is being "modulated" at some frequency.) Here the curvature is l/sin(cH), and hence 

X = /sin(ll)sin(cll)dH = .5;[cos((l-c)ll) - cos((l+c)H)]dH 
= (l/2(l-c))sin((l-c)H)-(l/2(l+c))sin((l+c)H) 
y = /cos(ll)sin(cII)dH = ,5/[sin((l +c)H) - sin((l-c)H)]dH 
= (1/2(1 -c))cos(( l-c)1 1) - ( 1 /2( 1 + c))cos((l + c)H) 

These are parametric equations for general hypocycloids and epicycloids ( [4], p. 312-314), with the 
substitutions t/2 = H and (l-c)/2 = b/a). When c>l, it is a hypocycloid; when c<l, an epicycloid; and when 
c = 1, a eycloid, as we showed in section 5. We can ignore the constants of integration - they just affect the 
location of die shape. 



As another example, take the Logo procedure ( [2], ch. 1) 

TO POLYSPI :S1DE 
FORWARD ;S1DE 
RIGHT 1 

POLYSPI SUM ;SIDE 1 
END 

which can also be written, following die discussion of section 7, and assuming a starting heading of 0, as 

rO POLYSPI 

FORWARD HEADING 

RIGHT 1 

POLYSPI 

END 

Hence K(H) = -1/H. Then from integration by parts: 

X = ; HsindDdll = sin(H) - Hcos(H) 
y = ; Hcos(H)dH = cos(l I) - Hsin(H) 

We cannot find an easy nonparametric forni in general, but note as H approaches infinity the second terms 

predominate, so 

X ~ -Hcos(H), y ~ -Hsin(H) 
r = Vlx" + y^] = H, 0 = arctan[tan(H)] = H 

and r = 0 is a polar expression for the curve. Thus as H becomes large die curve approaches the spiral of 
Archimedes ( (4], p. 317). We can show in general diat 
TO SPIRAL :N 

forward exponen t heading :N 
RIGH r 1 
SPIRAE :N 
END 

approaches as 1 1 becomes large, for N>0. 
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8.3. Connections to Fourier analysis 

The inverse of the curvature K is often called the "radius of curvature" R. Our equations for parameterized 

X and y can be written as 

X = -/ R(H)sin(H)dH 
y = -/ R(H)cos(H)dH 

and one can see these are just the formulas for the coefficients of sin(H) and cos(H) in a Fourier series 
expansion of R(H). 

There is another connection to Fourier analysis, however. Note if R(H) = Rj(H) + R 2 (H), then 

X = -/ R2(H)sin(H)dH - / R2(H)sin(H)dH 

y = -/ Rj(H)cos(H)dH - / R 2 (H)cos(H)dH 

and the resulting curve is the "vector addition" of the two original curves, as in DUOPOLY of ch. 3 of [2j. 
Now since we can approximate any curve sufficiently accurately with a sufficient number of terms from its 
Fourier series, we can approximate an arbitrary R(H) as the sum of weighted sines and cosines: 

R(H) = 2Japs(iH/c) + bjSin(iH/c)] 

where c is a constant representing how fine accuracy we desire in the approximation. Hence we can express x 
and y as 

X = -E./(axos(iH/c) + b.sin(iH/c)]sin(H)dH 
y = -E|/(a‘cos(iIl/c) + b|sin(iH/c)]cos(H)dH 

and the result is a weighted vector sum of hypocycloids and epicycloids, in the same way that MULTIPOLY 
in ch. 3 of [2] is the sum of circles of different sizes and curvatures. 

9. Conclusion 

We have shown a variety of methods for obtaining small-step turtle geometry programs from equations for 
analytic-geometry curves, and vice versa. Unfortunately, each metliod only works part of tlic time, and no 
guarantees can be given. See [5] for further student problems related to this material. 
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